SQUID Detection of Quantized Mechanical Motion 
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We predict that quantized mechanical motion can be detected by embedding a mechanical res- 
onator into a quantum SQUID. If the system is tuned to the regime when a plasma frequency of 
the SQUID matches the resonator frequency, the doubly-degenerate quantum level of the system is 
split by the coupling between the SQUID and the resonator. Observation of an avoided crossing as 
the function of external flux would be an unambiguous evidence of quantum nature of mechanical 
motion. We also investigate the conditions maximizing the level splitting. 
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The interest in nanoelectromechanical systems 
(NEMS) has been growing rapidly 0, Q because of their 
wide range of potential technological applications in 
detection and sensing, and their importance in testing 
fundamentals of quantum theory. The possibility to 
observe quantum mechanical motion of an oscillator has 
important implications in the understanding to which 
extent macroscopic objects obey the laws of quantum 
mechanics Q. NEMS can play also an important role in 
quantum computation where they have been proposed 
as qubits memory elements 0, @|, and quantum 
buses 0, 0|- Coupling of nanomechanical oscillators to 
a qubit has been also thoroughly investigated in the 
literature and many schemes of this kind have been 
proposed with Cooper-pair boxes d, 13. llH 
Josephson junctions (phase qubits) 
point contacts [lij], and quantum dots 




13, UJ, 

quantum 
16j . Recently the 
dispersive coupling of a NEMS to a Cooper-pair box has 
been realized [171 ]. 

Coupling of nanomechanical oscillators to SQUIDs has 
been recently intensively investigated 18, n| 2(J 21. 22. 
I23I ]. The high sensitivity of SQUIDs to tiny changes in 
the magnetic flux has suggested that the position of a 
nanomechanical resonator could be monitored by inte- 
grating the oscillator into the superconducting loop of 
a dc SQUID; indeed the transport properties of this su- 
perconducting circuit in presence of a uniform magnetic 
field depend on the position of the oscillator, since this 
position modifies the total area threaded by the flux. Re- 
cently this scheme has been demonstrated for the detec- 
tion of the thermal motion of a mechanical resonator in 
the classical regime [24| . Due to the high degree of control 
achieved on quantum SQUIDs, coupling nanomechanical 
resonators to SQUIDs is a promising scheme for observ- 
ing quantum effects in the motion of these oscillators as 
well. 

In this Article, we develop a protocol of detecting quan- 
tized mechanical oscillations with a SQUID coupled to a 
mechanical resonator. We show that the signature of 
this quantized motion is a splitting of an energy level 
associated with the SQUID due to the coupling to a res- 
onator. In the case of resonant coupling, this splitting 



can be detected by standard techniques developed for flux 
qubits 25]. We stress that achieving this regime is not 
straightforward as the typical frequencies of the SQUID 
(of the order of few GHz) and the oscillator (MHz range) 
do not naturally match to allow for a resonant behavior. 
In this Article we show that this is however possible by 
tuning the external magnetic field and the bias current 
with available experimental setups. 




Figure 1: A sketch of the device we study: one arm of a dc 
SQUID is free to oscillate in the plane of the SQUID itself. A 
uniform magnetic field B orthogonal to the SQUID plane is 
present and a dc bias current h flows through the device. The 
two Josephson junctions, whose phase drops are respectively 
71 and 72, are taken to be identical. 

The system we consider is schematically drawn in 
Fig. [TJ A dc SQUID is made of a superconducting loop 
of total area A with an arm of length I that can oscil- 
late freely in the plane of the loop. For simplicity we 
assume that only a single mode of oscillation with the 
frequency u) can be excited. This mode is described by 
the dynamical variable X representing the shift of the 
center-of-mass of the resonator (with the mass m) with 
respect to its rest position. The quantum effects related 
to mechanical motion of this oscillator appear at the scale 
of the amplitude of zero-point motion Xq = ^Jh/2muo. 
The SQUID further comprises two Josephson junctions 
of equal critical currents I c and shunting capacitances 
C; the typical energy scales related to the physics of the 
junctions are the Josephson energy Ej — hl c /2e and the 
charging energy E c = (2e) 2 /2C -C Ej, whose magnitude 
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depends on the geometry of the junctions. The typical 
time scale for the dynamics of the junctions is set by the 
inverse plasma frequency lu p \ = y/2EjE c /h. The dynam- 
ics of the SQUID is described by the two gauge-invariant 
phase drops 71 and 72 across the two junctions. 

The coupling of the SQUID dynamics to mechanical 
motion is provided by the position dependence of the 
magnetic flux threading the SQUID loop. The two phases 
are constrained by the requirement that the supercon- 
ducting order parameter is single valued, 



7l — 72 = 27T 



(1) 



where $0 = h/2e is the flux quantum; the total flux $ is 
the sum of three contributions. The first one comes from 
the external bias <E> e = BA = ^o^e, while the second 
one <& m = BIX depends on the position of the mechani- 
cal resonator and provides the coupling between the me- 
chanical resonator and the SQUID. If the circuit also has 
non-negligible self inductance L, the third contribution 
to the total flux <E> comes from the current circulating in 
the loop; if I\ and I2 are the currents flowing through 
each junction, the self-induced flux reads L(I\ — /2V2. 
This device has three degrees of freedom; if L = 0, then 
the constraint expressed by Eq. ([T]) reduces the number 
of degrees of freedom from three to two. In the following, 
we assume that the dissipation effects arc negligible. 

Quantum coherence in the motion of the mechanical 
resonator can be detected via spectroscopic measure- 
ments on the quantum SQUID. Indeed, the energy lev- 
els associated with the SQUID degrees of freedom are 
shifted due to the coupling to the oscillator. However, 
this shift is of the second order in the coupling and for re- 
alistic experimental parameters is very small. We employ 
therefore below a different scheme, which provides the re- 
sult which is of the first order in the coupling. One first 
chooses an appropriate set of values for the externally 
controllable quantities if, and $ e such that the system 
can be trapped in a minimum in which one of the eigen- 
frequencies of the electromagnetic modes is the same as 
the frequency of mechanical oscillations (to be referred 
below as degeneracy condition); then one moves slightly 
away from this degeneracy condition by changing the re- 
maining external parameter. A plot of the energy levels 
as functions of the external parameters should display an 
avoided level crossing, which is a clear indication that the 
SQUID is coupled to a coherent quantum system. 

It is useful to describe the system in terms of three 
dimensionless variables 7 = (71 + 7a)/2, <fi — $/$o and 
£ = BIX/Qq. The potential energy of the system reads 



U = Ei 
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where we have introduced the parameters 
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The parameter si is proportional to the flux threading 
the area swept by the mechanical resonator and plays the 
role of a coupling parameter, as we show below. Typ- 
ical values of si — 4 x 10 -5 can be obtained assum- 
ing I c — 1/iA, lu = 1GHz, / = lfim, Xq = lOfm and 
B = 0.1T. When the temperature is low enough for the 
system to reach the quantum regime, the coordinates typ- 
ically oscillate around a minimum of the potential corre- 
sponding to the values 7, <j) and £. We can approximate 
the dynamics as a three-dimensional harmonic oscillator, 
with the energy being a quadratic form, 
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where the coordinates qi read 



<h = 7 - 7; 92 = n 

and 
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(We have introduced r = cos7cos(7T0), s = sin7sin(7r0) 
and 17 = uj/oj-px). The coupling V23 between the me- 
chanical resonator and the SQUID is proportional to the 
ratio si / (3l (for devices with a low self-inductance this 
behavior no longer holds, see below for discussion); the 
decoupled regime can be obtained for either B — ► or 
L — > 00. The coordinate Q\ =7 — 7, corresponding 
to one of the electromagnetic degrees of freedom, oscil- 
lates with the frequency uj p \ y cos 7 cos(-7r</>) ; therefore if 
a minimum (7, <f>, £) is such that 



COS7COs(7T</>) = fl 2 , 



(7) 



the frequencies associated with the motion of average 
phase drop 7 and oscillator motion coincide up to correc- 
tion of second order in si; we show below that the equal- 
ity is indeed exact to all orders in the coupling. Since in 
the low-inductance limit the coordinates <p and £ are not 
independent, it is convenient to switch to the basis where 
the submatrix corresponding to £ and 4> is diagonalized. 
The analytic expression for the matrix in the new basis is 
quite cumbersome; below we give this expression in the 
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limit in which the degeneracy condition {7J holds, 
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with C — \ / 2tt 2 £/ 2 9 2 + 1. The first and second row and 
column correspond to the phase drop 7 and to the me- 
chanical degree of freedom, respectively. The parameters 
lb and <& e corresponding to the degeneracy condition are 
found if one solves the equation set, 



h/Ic — 2sin7cos(7r0) 

2/3^ 1 (<j> - cj) e - £ ) = - cos 7 sin(7r<; 



£ = 47r^/£ A (0-&-£- 
cos7cos(7r0) = 9 2 , 



(J3 L I b )/(4I c )) 



(9) 



where the first three lines are satisfied by stationary 
points of the potential energy ^ and the fourth line is 
the degeneracy condition. In a minimum one must have 
V > 0. The unknowns are the coordinates (7, of the 
minimum and <fr P ; we use lb to tune the system to the 
degeneracy point [261] . 

Next, we quantize the system. The Hamiltonian reads 

3 1 

! i& 4 V W i W 3 

_ (10) 
with LUi — ^ v \\fV[i- When the condition ([7]) is ful- 
filled, the first excited levels 1 100) and 1 010) are quasi- 
degenerate and the Hamiltonian (fTU|) restricted to their 
subspace becomes 
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H = 
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A = ^ pl ^' 7 . (11) 



Fig. [2] shows the dependence of the dimensionless level 
splitting \/fiuj v \ on the dimensionless external flux <f> e . 
The maximum level splitting is obtained for a specific 
value of the external flux and depends very weakly on 
the loop self-inductance; however the value of the flux 
at which the maximum is achieved, does; see below for 
further discussion. The magnitude of the splitting is of 
the order of 10 _5 ?ia;pi. The value of the ratio 9 — lo/lo p \ 
plays an important role: lower values correspond to big- 
ger maximum splittings and are thus preferable. 

To enable the detection, the potential well formed at 
the chosen minimum must be capable of containing quan- 
tum states. The number of bound states can be esti- 
mated by the ratio between the energy difference AU 
between the minimum and the closest saddle point and 
the energy level separation, which is roughly Too, 
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Figure 2: Dimensionless splitting \/Tiuj p \ versus the external 
flux (f) e for the coupling parameter stf — 4 x 10~ 5 . The self- 
inductance parameter /3l is 10~ 4 for the left panel and 1 for 
the right panel. Dotted, dashed, and solid lines correspond 
to f2 = cj/ojpi = 0.1,0.5 and 0.9, respectively. The plots are 
periodic, the period being one flux quantum. 
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Figure 3: (a) the dimensionless depth Au of the minimum 
featuring maximum splitting A is plotted as a function of the 
ratio SI = w/w p i, for different values of the self-inductance 
parameter: /3l = 10 -4 (solid line), /3l = 1 (dashed line) 
and (3l = 10 (dotted line), (b) the maximum value of the 
dimensionless splitting \/hw p i is plotted as a function of Q; 
here = 1, but the result is independent of /3l- 



Fig. [3^ displays plots of the quantity Am, defined by 
Eq. (fl"2|) . as a function of 91 for different values of the 
sclf-inductance parameter 0l\ one sees that small values 
of 9 correspond to very shallow minima in the poten- 
tial energy, whereas the deepest minima correspond to 
devices in which lu = lu p i. However Fig. [3Jd shows that 
for these devices no gap is expected in the first order in 
the coupling; this rather surprising fact can be further 
illustrated by considering the dependence of the matrix 
element VL on s — sin7sin(7T0) in Eq. ©: for 0=1, 
one has cos7cos(7r0) = and thus s = = A. Because 
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of this trade-off, the optimal condition corresponds to an 
intermediate case. Note however that Am depends on 
Ej and E c only via the ratios Q — lo/lo v \ and Ib/I c and 
thus can be tuned independently of Ej/E c . Therefore a 
bigger minimum depth can be obtained by designing the 
SQUID so that the ratio Ej/E c is large. 

We now comment on the role of the self-inductance. 
The curve in Fig. [3b is very little affected by the value 
of Pl, implying that this parameter is not relevant for 
improving the splitting A; on the contrary, the curves 
in Fig. [3^ show that the dimensionless depth Au of the 
minimum well is drastically reduced by increasing the 
self-inductance of the loop. Thus loops of smaller self- 
inductance are preferable. 




Figure 4: Plots of the solutions of Eqs. giving maxi- 
mum value of the dimensionless splitting A/?kj p i; left panel 
(Ib/Ic, 4>e, fi), right panel (</>, 7,f2). A projection of the points 
onto to the horizontal plane is also displayed for reader's con- 
venience. 

Finally, Fig. [4] represents numerical solutions of 
Eqs. used in the previous figures. The right panel 
shows the values of the dimensionless external parame- 
ters lb / I c and <p e that corresponds to the maximum value 
of the gap A for values of f2 between and 1, while the 
left panel shows a plot of the coordinates (7, 4>) of the 
best minimum. As expected, the results are symmetric 
with respect to simultaneous current inversion If, —> —lb 
and flux reflection 4> e — > I — <j) e ; this can be seen from 
Eqs. ©, where in this case if (7, 4>, £, Ib/I c ) is a solution, 
then also (n — 7, 1 — 0, — f, —Ib/I c ) solves the equations. 

In conclusion, we found that quantum oscillations of a 
NEMS embedded into a SQUID can be detected by spec- 
troscopic measurements in the regime when one of the 
plasma frequencies of the SQUID matches the frequency 
of the mechanical resonator. This frequency matching is 
possible with the current experimental techniques, and 



the scheme has two parameters (external flux and exter- 
nal current) to simultaneously tune the system to the 
vicinity of the degeneracy point and perform spectro- 
scopic measurements around this point. Measurements 
of splitting of the degenerate doublet state displaying an 
avoided crossing would be an unambiguous evidence of 
quantum nature of mechanical vibrations. 
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